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1 Introduction 


The quantum graph stated is a pure quantum n-qubit state whose joint stabilizer can be 
written using a symmetric matrix whose elements are Pauli matrices. The Pauli matrices 


are: 




Z 



Y = iXZ, 


where i = \/^. More precisely, in [T], graph states in n qubits are defined as the unique, 
common eigenvector in (C^)” to the set of independent commuting observables: 


Ka = XaZM^ =Xa\[Z^ 
b&Na 

where Aa = I ® I ® ® A® I ® ® 1 1 Na are the neighbours of a in the associated graph 

(see below), and where the eigenvalues to the observables Ka are +1 for all a = 0,..., n — 1. 
A natural generalisation is given by 

Xq (JaZj^f^ (T a 11 Zi) 

beJVa 


where a G {X,Y}. For the mixed graph states of this paper we require this extra generality. 

We represent a pure quantum graph state by a simple graph, hence the name graph 
state. For instance, for n = 3, the simple graph with vertices V = {0,1,2} and edges 

0 1 1 \ 

10 0 and can be used to represent the 

1 0 0 y 

pure quantum state of 3 qubits that is jointly stabilized (with eigenvalue 1) hy X ® Z ® Z, 
Z®X®I, and Z®I®X. This set of n = 3 joint stabilizers can be written in matrix form as 
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E = {01, 02} has an adjacency matrix A = 


1 



/lll\ ! X I I \ ! I Z Z \ ! X z z \ 

A = / XI * Z / I j = I Z X I , 

yiiiy \i I X j \z I I j \ z I X j 

where A is the element-wise product, is the 3 x 3 identity, and the connection with the 

0 1 1 \ 

1 1 1 , we get 

1 1 l) 

/lll\ / X I I \ / I Z Z \ / X z z \ 

A = I x l XI y^lzZz\ = lzYZy 

\l I i J \l I X j \Z Z Z j \ Z Z Y j 

As the system is stabilized by any of the n matrix rows, it is therefore also stabilized 
by any product of the rows. In general, this set of stabilizers forms a stabilizer code with 

0 1 1 \ 

I 0 0 , we have that 

I 0 0 y 

{X ® Z ® Z){Z ® X ® I) = (Y <Z)Y Z) I) is also a joint stabilizer of our graph state, as is 
{ZZIZX){XZZZZ){ZZXZI) = -{X ZYZY). 

Given the matrix A of stabilizers the implicit assumption above is that these elements 
are stabilizing a pure quantum state, i.e. that the n qubits are not entangled with an 
environment. Given this assumption then the rows of A must commute, otherwise they 
cannot dehne (stabilize) an n-qubit pure state. Let l"^) be this pure state and let U be any 
row of A. Then U \'ip) = \'ip), where \'ip) is unique, the density matrix of the graph state is 
p = 1'^) ('^1, and UpU^ = p. 


2 ” elements, i.e. a stabilizer group, S. For instance, for A = 


adjacency matrix, A, should be clear. As another example, for A = 


Lemma 1 Let p = \%Ij) {%Ij\ he the density matrix of a pure graph state, and U as above. 
Then, as U is taken over all rows of A, then UpU^ = p determines p uniquely up to local 
Pauli unitary equivalence. 

Proof. Since G is a tensor product of local Pauli matrices, U If)) {iplW = It/’) (t/’l implies 
G It/’) = ± It/’)- This means that It/’) is a common eigenvector of the set of matrices, U, with 
eigenvalue ± 1 . Let P be an n-fold tensor product of Pauli matrices. Then it is straightfor¬ 
ward that UPpP^U^ = PUpU^P^ = PpP\ so It/’) is dehned by U up to local Pauli unitary 
equivalence. □ 

The symmetric matrix A is completely characterised, up to, but not including, commu¬ 
tation, by Ai = A + uiln, where Aa is an additive matrix over F4, is the n x n identity 
matrix, a; is a primitive generator of F4, i.e. = 1 , and where / —)■ 0 , Z —)■ 1 , X —)■ ta. 



j 0 

1 

1 \ 



1 

1 \ 


y —)■ For instance, for A = 

1 

1 

1 

) *- 4-4 — 1 

1 

0.2 1 

. The set of 2 ” sta- 


\ 1 

1 

1 / 


\ 1 

1 

0.2 / 



bilizers then become the set of 2” codewords of an F4 additive code (where multiplication 
of rows of A becomes F4 addition of rows of Aa). This code is self-dual with respect to 
the Hermitian inner product because the matrix is symmetric with u or on the diagonal. 
Observe, however, that this map from stabilizers to F4 does not take into account matrix 


2 



commutation as it does not capture the minus sign of, for instance, XZ = —ZX^ as this 
translates to tc.l = l.tc. 


In this paper we define and characterise mixed graph states, by considering square ma¬ 
trices A that are not, in general, symmetric. In previous work [ 3 ] the resulting F4-additive 
codes (not, in general, self-dual) have been classihed and referred to as directed graph codes. 


For instance, the adjacency matrix A 


0 1 0 \ 

0 0 1 is not symmetric and describes a mixed 
0 10 / 


graph comprising an undirected edge 12, and a directed edge 0 —1. 
As before we can map such a matrix to a matrix of stabilizers, thus: 


Z X ) 

But what is a natural quantum interpretation of such a matrix? Rows 0 and 2 , and 1 
and 2 commute, but rows 0 and 1 anti-commute. A quantum interpretation is only possible 
if all rows pairwise commute, so we choose to extend the rows by extra columns until the 
rows pairwise commute. Implicit in this column extension choice is the assumption that we 
add environmental qubits as opposed to, more generally, qudits or continuous variables. 

In this paper we show how such matrices can dehne a class of mixed states, which we 
call mixed graph states. The products of the rows of the matrix still form a group, S, up to 
a global constant, but now the group does not have a joint eigenvector, i.e. we don’t have a 
pure state, since some of the rows anti-commute. 

Due to anti-commutativity, some of the elements of the group are only dehned up to a 
global constant of ±1 (there is no natural ordering on the rows of A). For instance, for the 
matrix A associated with our example above, [X Z) Z Z) I){I Z) X ® Z) = —{IZXZ) Z){X Z) 
Z (g) /), as the two rows anti-commute. So S contains ±z(X g) R Zi Z). 

We can than dehne a density matrix, p, which is stabilized by all members of S. Specih- 
cally, 

sps^ = p 'is & S. 

More generally, p is stabilized by any as, where |q:| = 1 , which is why we don’t have to 
concern ourselves with global phase constants in members of S. 

Throughout this paper we dehne the pure state vectors using a Boolean function no¬ 
tation [ 5 ]. Specihcally, we associate with our n-qubit graph state an n-variable homoge¬ 
neous generalised quadratic Boolean function p : F2 —)■ Z4, where p{xq,xi, ... ,Xn-i) = 
Ylj<k‘^^jkXjXk + Aj. jXj, and A is the modihed adjacency matrix of our graph state with 

mixed graph is a graph where some of the edges may be directed. In contrast, all edges of a directed 
graph are directed. 
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j ^ jk 

elements defined as Aj^ = \ 0 ,j = k and aj = X , with T the adjacency matrix of 

1 ,j = k and aj = Y 

the graph. Then one can write the graph state \^jJ) as \^jJ) = 

In later sections, we refer to the nodes where aj = X as white nodes, and to the nodes 

where (Xj = X as red nodes. 


Example 1 We show here the mixed graph state assoeiated with the directed triangle: 



2 


The stabilizer group S is generated, up to ±1 constants, by: 

X ® Z 
I ®x® z 
z ®x 


Because some of the rows of A anti-commute, instead of commute, there is no pure quantum 
state 1'^) such that (X ® Z 0 /) |'^) = (/ ® X ® Z) |'^) = (Z ® ® X) |'^) = |'^). 

Instead we can associate S with the mixed graph state whose density matrix is: 


Po = 


/ 1 
0 
0 

—i 

1 

0 

0 

V ^ 


0 

1 
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0 

0 

-1 

—i 
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0 

i 

1 

0 

0 

—i 

1 
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i 

0 

0 

1 
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0 
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1 

0 

0 

—i 

1 

0 
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0 

-1 

i 

0 

0 

1 

i 
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0 

i 

1 

0 

0 

—i 

1 
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\ 

0 
0 

-1 
—i 
0 
0 
1 


We will discuss in later sections how we obtain this density matrix. It is easy to check that 
{X ® z (Z I)pq{X ® z ® = {I ® X ® Z)pq{I ® X ® Z)"^ = {Z ® I ® X)pq{Z ® I ® X)"^ = Po. 


Po is not the unique mixed graph state associated with S, however: Other density matrices 
Pj such that 

{XZ)ZZl)pjiXzZZlY = {I ® X ® Z)pj{I ® X ® Z)"^ = {Z ® I ® X)pj{Z ® I ® X)"^ = Pj 
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are: 
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—i 
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0 

0 

-1 

1 

—i 

—i 

0 

0 
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0 

i 

0 
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1 

0 

-1 
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,P2 = 
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-1 
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0 
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0 
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0 
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-1 


V 

i 

0 

i 

0 

0 

-1 

0 

1 




V 

i 

i 

0 

0 

0 

0 

-1 
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There are also another 3 matrices which are locally equivalent to these 3 (see lemma 24). 
Thus we dehne the mixed graph state associated to the graph as 
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P ~ 

i =0 

where Cj > 0, — 1) i-®- convex sum of these density matrices. Note that, 

although only at most 3 of them are locally inequivalent, we include them all in the general 
sum p, since changing pj for a locally equivalent p'- does not necessarily give local equivalence 
in p. 

Each of the pj can be associated with a pure graph state in 4 qubits. For instance 


Po 


\ '^j2{xoX2+X\X2)+X2^ ^j2{xQX2+X\X2)-\-X2 



j2(xoX2+X-iX2+X-i+X2)+X2 'j l^j2{xQX2+XiX2+Xi+X2)+X2 


is the density matrix resulting from a measurement of qubit 3 of the pure graph state 
j 2 {xox 2 +xiX 2 +xiXi+x 2 xP)+x 2 terms of Boolean functions, this would correspond to the evalua¬ 
tion of X 3 at 0 and 1). Similarly, pi and p 2 are the density matrices resulting from a measure¬ 
ment of qubit 3 of the pure graph states i‘^{^oxi+xox3+x2xz+xo)+x2 j 2 (xoxx+xix 2 ^xoxz+xxxs)+xi ^ 

while the remaining pj will be given by the addition of linear terms to the polynomial. 


2 Mixed graph states 

/ X z I 

To recap, a matrix such as ^ = j I X Z 

\ I Z X 

muting, but can still be interpreted as a quantum object by making it part of a larger fully 
commuting matrix, i.e. where we choose the environment appropriately, and this will imply 
that our quantum object is a mixed state, in fact a mixed graph state (in this paper, we limit 
our study to the addition of environmental qubits as opposed to, more generally, qudits). 

Let G be the mixed graph dehned by adjacency matrix A, and the undirected graph 
that has adjacency matrix T = A -|- Thus Gfe is the simple graph obtained from G by 
erasing all undirected edges of G and making all directed edges of G undirected. 


^ has rows that are not fully pairwise com- 
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What is the minimum number of columns, e, necessary to append to A so as to make 
all rows pairwise commute? We can then add the same number, e, of rows to the matrix 
so as to make it square again, making sure that all resultant rows pairwise commute with 
the previous rows and with each other. For instance, we can extend our example by e = 1 


columns and rows to get A1 = 


( 

Z 

\ 


f ^ 

z 

I 

^ \ 

A 

X 



I 

X 

z 

X 


I 



I 

z 

X 

I 

\ Z I I 

X 

J 


\ Z 

I 

I 

X J 


which is fully 


commuting. So for this example the minimum number of columns and rows necessary to 
make the matrix fully commuting was e = 1. Finally, as this matrix is fully commuting we 

/ X Z I Z \ 

can, by suitable row multiplications, recover its graph form = “f "5 ^ { 

1 Zj Ji. 1 

\ Z I I X J 

We consider the pure graph state described by A^ to be a parent graph state, being a 
parent of the mixed graph state described by A. We then denote the resultant mixed state 
given by tracing out the environmental qubits as its child. 

! X Z I \ /oio\ 

For our example, ^ = I / X Z , we have A = I 0 0 1 , and the minimum 

\i z X j yoloy 

number of columns and rows necessary to add to A to make it fully commuting is 


e = 


1 

2 


rank(F) 


-rank 

2 


0 10 
1 0 0 
0 0 0 


1 


In general: 


Lemma 2 (special case of |i4j) The minimum number, e, of columns and rows required 
to be added to A to make its rows pairwise fully commuting is given by 


e = 


1 

2 


rankfr). 


( 1 ) 


We refer to e as the mixed rank. 


For our example we extended to a 4-qubit pure graph parent state where the 4th qubit 
is part of the environment. This pure graph state can be written using Boolean function 
notation as 




Tracing out the 4th qubit, means summing the projectors obtained by fixing X 3 = 0 and 1, 
respectively. We obtain the projectors |0o)(0o| and |0i)(0i|, where |0o) = W(^—iyoxi+xiX 2 ^ 
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|0l) = ^(^ — l'^^0Xl+XiX2+X0^ 


p 


I0o)(0o| + |0l)(01 


/ 


-1 -1 

1 1 

1 1 

-1 -1 


1 

0 

1 

0 

-1 

0 

1 

0 

1 


1 

1 

0 

1 

0 

1 

0 

-1 

0 


-1 

-1 

-1 

1 

-1 

-1 

1 

-1 

0 

-1 

0 

1 

0 

-1 

0 

-1 


1 1 

1 1 

1 1 

-1 -1 

1 1 

1 1 

-1 -1 

1 1 

1 0 

0 1 

1 0 

0 -1 

1 0 

0 1 

-1 0 

0 1 




So we can interpret the stabilizer matrix A as representing the mixed qnantnm state p. 
However this is not a complete interpretation for two reasons: 


• One can extend by more than one colnmn/row and still obtain a fnlly commnting 
matrix. 


• The matrix A can be extended in mnltiple ways, leading to mnltiple, possibly ineqniv- 
alent, mixed states p. 


We shall, initially, avoid the hrst issne by stipnlating that we only extend by the minimnm 
possible nnmber, e, of colnmns/rows, as given by Q. This is compatible with the pure graph 
state formulation as e = 0 forces the state to be pure. As an example of the second issue, ob¬ 


serve that it is equally valid to extend ^ to = 


from which, by multiplicative row operations, we obtain A^ = 


/ 

X 

\ 


f ^ 

z 

I 


A 

z 



I 

X 

z 

z 


I 



I 

z 

X 

I 

\ I Z I 

X 

J 


\I 

z 

I 

X J 


I 

I 

V I 


I 

X 

z 

z 


I 

z 

X 

I 


I \ 

z 

I 

X J 


It is 


clear (since the former is a connected graph and the latter is not) that the parent graph 
state described by this A^ is locally inequivalent to the previous parent. Therefore the resul¬ 
tant mixed state, p = |0o)(0o| + obtained by tracing out the 4th qubit, is different 

from the previous mixed state. In general we shall have multiple parent graph states leading 
to multiple density matrices, so our overall density matrix shall be a convex sum of these 
density matrices. 


Dehne the codespace (i.e. the stabilizer group S) of the n x n matrix. A, to be the 2"' 
stabilizers formed by products of one or more of the rows of ^ - remember that, as some 
of the rows of A anti-commute, some of the members of S are only dehned up to a global 
multiplicative constant of ±1. Likewise, the codespace of A'^ comprises stabilizers, but 
now all rows of A'^ commute so the global constant is always 1. 
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as 


We refer to H = ^ ^ ^ Hadamard matrix, and to 

the negaHadamard matrix. Both are unitary. 

Lemma 3 The Pauli group, {I , X, Z,Y} , is conjugated (up to a factor of ±1) by the group 
generated by Specifically, 



I 

H 

N 

iV2 

NH 

HN 

X 

X 

Z 

-Y 

z 

X 

Y 

z 

Z 

X 

X 

-Y 

-Y 

Z 

Y 

Y 

-Y 

-Z 

-X 

Z 

-X 



For instance, X = HZW, —Z = NYN\ and —X = N‘^Y{N‘^y. Lemma implies that, by 
conjugation, to within a factor of ±1, we can permute the elements of the stabilizer matrix, 
column-wise. There are 3! = 6 such permutations. 

The following is well known (proof omitted): 

Lemma 4 Consider a guantum system comprising a local system, L, possibly entangled 
with an environment, E, and let pi be the density matrix that defines the local system. Then 
unitary conjugation, Ue{L x E)Uy, on the environment leaves pi unchanged. 


Lemma 5 The extension columns and rows do not depend on the direction of the arrows, 
although the parent graphs depend on the direction of the arrows. 

Proof. The extension columns and rows depend exclusively on the commutativity or 
anti-commutativity of the stabiliser basis, which is not dependent on arrow direction. □ 

3 General Density Matrix Stabilized by A 

A density matrix of n qubits can be written in a Pauli basis, i.e. where the basis elements 
are tensor products of elements from {/, X,Y, Z}, i.e. 

P = where k = {ko,ki,... ,kn-i) e Z2,(Tk = (Tkj, 

do = /, di = X, d 2 = F, ds = Z, and Xlfcezj WA <1; ^ 

where the Ok are real. Necessary and sufficient conditions for p to be a density matrix, i.e. 
derived from a statistical sum of pure states by tracing out the environment, are: 

• p is Hermitian (i.e. equal to its transpose conjugate). 


• p is positive semi-dehnite. 

• Tr(p) = 1. 










We refer to P = {a^ | k G ZJ}, as the set of Pauli codewords (stabilizers) of size 
|P| = A"'. Let A be the nx n matrix for a mixed graph state. Then A"^ is the matrix for the 
mixed graph state obtained from the graph for A by reversing all arrows. As before, we say 
that A generates a Pauli subgroup, S' C P. We will now show that A'^ generates the dual 
group of S, namely S'-*- C P. 

Theorem 1 The set of length-n Pauli words (stabilizers) that commute with all rows of A 
forms a multiplicative group, S-*-, of size 2”, being precisely the words generated by A^, to 
within ±1 constants. 

Proof. It is straightforward to show that every row of A'^ commutes with every row of A. 
To show that there are no more, observe that each row of A^ can be viewed as a restriction. 
Only 4"/2^ (i.e. half) of all possible length-n Pauli words commute with a given row. The 
rows of A^ are independent so the n rows of jointly commute with only 4"'/2" = 2"' Pauli 
words. But the rows of Af’" generate 2" words, which must be distinct and, therefore, must 
be all of them. □ 

See section [3 for a discussion of the ±1 coefficients. 


Graph states can be also described as the sum of its stabilisers, as seen in [?]. In order 
to state a similar result for mixed graph state, we need the following lemmas: 

By the support of a matrix, M, we mean the set of non-zero element positions of M, i.e. 
support(M) = {{i,j)\Mij ^0}. 

Lemma 6 The elements of the stabilizer group, S'-*-, of a mixed graph have non-intersecting 
support. 

Proof. By inspection X and Y have the same support, as do I and Z. But the support 
of X and Y is non-intersecting with that of I and Z. Moreover if two matrices dk and dfc/ 
are both tensor products of Pauli matrices, and if they differ in at least one tensor position, 
where one matrix has X or Y, and the other matrix has I or Z, then dfc and dk' must also 
have non-intersecting support. The matrix A only has X and/or Y on the diagonal and I 
and/or Z off the diagonal. Let A have rows numbered 0 to n — 1. Let R and R' be two subsets 
of {0,1,..., n — 1} and consider the matrices, dk^ and , being the product of the rows 
of A indexed by R and R', respectively. Then dk^ and ak^, have X or F at tensor positions 
indexed by R and R', respectively, and J or Z at all other positions. So, unless R and R' 
are equal, they must be non-intersecting in at least one tensor position, so therefore and 
have non-intersecting support. The elements of the stabilizer group of the mixed graph 
are obtained by ranging R over all subsets, i.e. over all stabilizer codewords. Therefore any 
two of them must be non-intersecting. □ 

For a given stabilizer group S', the condition sps^ = p Vs G S' is equivalent to the condition 
sp = ps \/s ^ S'. The Pauli basis description of ([^ imposes the condition that p = ^kd'k, 
where the dk are tensor products of Pauli matrices, such that sdk = dkS Ws E S. In other 
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words, p is the sum of tensors of Pauli matrices that commute with all s G S. By theorem 
the Pauli subgroup of matrices that commute with all elements in the stabilizer group S 
is the stabilizer group, S-^, of the mixed graph with the arrows reversed. The intersection 
of S and S'-*- is then the row space of A associated with the vertices of the mixed graph, G, 
that are isolated in Gb. But what choices do we have for the angle of in ([^, such that p 
is also Hermitian (as a density matrix must be)? Let a*, = 0 ,^ 0 ,^ such that |afc| = 1 and 
is real. Lemma implies that each component ak^k of P must, itself, be Hermitian, so this 


hxes hk precisely, to within ±1. For example, consider A^ 


X Z 
I X 


Then contains 


(712 = ±i(x (g) y) = ± 


( 0 
0 
0 

V-1 


0 

0 

-1 

0 


1 \ 
0 
0 
0 / 


which is not Hermitian, so we choose 012 = i or 


ai 2 = —i in this case. In the sequel we deal with subgroups, S and S-^, of the tensored Pauli 
group, of size 2 "' (to within ±1 coefficients). So instead of indexing by A; G Z 4 , we index 
using K C = {0,1,..., n — 1} relative to the rows of the group generating matrices A 
and A^^ respectively. Specifically, 


P=Y1 ^ ^kAk, 

fcezj kcz „ 

where the bx G {±l,±i}, and A]^ = Ylh&K^li^ where A]^ is the hth row of A'^, and 
= -AJX if j' > j, i.e. we fix the product ordering with lowest indices on 
the right. Moreover, for V obtained from A^, bx is unique if the elements in K are pair¬ 
wise commutative, and is otherwise only defined up to ±1. For instance, ii K = {2,4,5} 
then and bx is unique only if the elements in {2,4,5} are pairwise 

commutative. 


4 The maximum commutative subgroups of 

We now show that any commutative subgroup of the group of operators generated by the 
rows of A^ = JA is contained in a commutative subgroup of maximum size 2”“®. Firstly, we 
show how the commutativity of the Pauli operators of a mixed graph state can be expressed 
in terms of the associated simple undirected graph, Gb, and its corresponding adjacency 
matrix, F = Fg^. Then, in section we show that the density matrix of any child of a 
pure graph state parent can be represented as a weighted sum of the terms of a maximal 
commutative subgroup of S'-*-. 

Lemma 7 Let Af he row j of A^. Let J,K C {0,1,..., n — 1}. Then Af commutes with 
Ax iff the number of elements of K Ci Mj is even. Furthermore, Aj commutes with Aj^ iff 
the sum over J of the number of elements of K r\ Mj is even. 
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Proof. Af anticommutes with any Aj^ such that k & K r\ N'j, that is, any k in the 
neighbourhood of j, thus AfAj^ = —Aj^Af. If there are an even number of neighbours 
then the minus signs cancel each other out so that Af commutes with .4.;^, otherwise Af 
anticommutes with A^. More generally, Aj commutes with A^ iff, after passing, for all 
i G J, Aj through an even number of minus signs are generated, otherwise 
anticommutes with A^. □ 

Example 2 Let Gb be the simple graph 01,04,12,23,34,14- Then A^ commutes with 
(since they are independent, so the intersection is empty), and both commute with Aj^ = 
*^{12 4 } ~ A):A 2 Af, since 0 has 1 and 4 os neighbours, and 3 has 2 and 4 os neigh¬ 
bours, so both have two elements in the intersection with K = {1,2,4}. Note that e.g. 
*^0 (*^ 4 Mg = —A^A^AqAi = +A^A 2 AiAq. We also have that 

A^A^ commutes with A^AiAq, since both 2 and 3 have one neighbour in K = {0,1,4}, 
so the sum over J = {2, 3} is even. 


Let Tkj be the {k,j)th element of L = AA^. 


Lemma 8 Af commutes with Af iff ®keK^k,j = 0, where '©' is the binary sum. Further¬ 
more, Aj commutes with Af iff 0jgj = 0. Let vk = ("Oo, 'Oi,..., Vn-i) G be 

such that Vk = I if k E K and Vk = 0 otherwise. Then the property that Aj commutes with 
Af translates to the condition = 0. 

Proof. We hrst prove that Af commutes with Af iff ~ lemma 0© 

commutes with Af ^ the number of elements of K flA/}- is even, i.e. there are an even number 
of rows k E K with Tkj = 1, i.e. Tkj = 0. Now we prove that Af commutes with Af iff 
0ieJ 0fcei^ Tfcj = 0: By lemmapT 


even, i.e. 


Af commutes with Af iff the sum over J of the number of 
there are an even number of rows k E K such that L kj = 1, i.e. 


elements of K OAfj is 

= 0. Finally, note that vkT = {^k&K^ 
from which it follows that VK^vf = 0ie J 0 keK ^k,j- 

As r is symmetric, we get that uxTuf = vjFvf = 0 , which ensures that commutativity 
is a symmetric relationship. 


k&K • • • ) ®keK rfc,n-l), 

□ 


Example 3 For T = 


( 0 
1 
0 
0 

VI 


1 \ 
1 
0 
1 
0 / 


Mg commutes with Af since Fg^g = 0, and both 


Af and Af commute with AfAfAf, since 0{i_2,4} Ffc,o = l©0©l = 0 and 0{i^2,4} = 

0 © 1©1 = 0 . We also have that AfAf commutes with AfAfAf, since 0 jg {2 3 } 0fce{o 14 } ^^4 
(0 © 1 © 0) © (0 © 0 © 1) = 0, i.e. vkFvj = 0, for vk = (1,1, 0, 0,1) and vj = (0, 0,1,1, 0). 


If Af commutes with Af and Af, commutes with Af, and J fl J' = 0, then AfAf, 


commutes with Af. Similarly, if Af anti-commutes with Af and Af, anti-commutes with 
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and Jn J' = 0, then A^Ajji commutes with Aj;^. If one of them commutes and the other 
anti-commutes, we get anti-commutativity. This implies that subgroups can be combined to 
generate other subgroups in an easy way. 

Lemma 9 The subgroup structure of S~^ (and therefore also of S) is independent of the 
direction of the arrows in G: those graphs, G, sharing the same Gb will have the same 
commutative subgroup structure (although the precise group members differ). 


Proof. The subgroup structure is only dependent on Gb, which is the same regardless of 
arrow direction. □ 


From lemma 19 
members of S 


any two vectors v, v' G F 2 that satisfy vTv''^ = 0 dehne two commuting 
From here it follows that, if F has rank n — t, then comprises 2* copies of 


a group that is isomorphic to S^, where S'-*- is generated from the rows of A^, where A^ is 
obtained from A^ by deleting t dependent rows and corresponding columns of A~^. We then 
consider the maximum commuting subgroups of S'-*- instead of S'"*". Equivalently we reduce 
F to F, being a.n — txn — t matrix of maximum rank, and look for two vectors v,v' E 
that satisfy vTv'"’" = 0. More generally, by linearity, there exist size 2^ subgroups of 
(linear subspaces). Lb, as generated hj h x n binary matrices, B, that satisfy BTB^ = 0. 
Each such linear subspace. Lb, dehnes a commuting subgroup of S'-*-. 


Lemma 10 Let F have rank n — t, and let T be an — txn — t maximum rank matrix 
obtained from F by removing t linearly dependent rows and corresponding columns. There 
are multiple choices for T. Then each commuting subgroup, P of S^, as described by some 
B satisfying BTB^ = 0, is in one-to-one correspondence with a commuting subgroup, P, of 
S'-*-, as described by some B satisfying BTB^ = 0, where \P\ = 2*|P|. 


Proof. Let Q be a set of t linearly dependent rows of F. For some u, u' G F 2 let = 0. 

Then u,u' G for some 13. Then 3w,w', uniquely, such that uT = wT and u'T = w'T, 
where w and w' are both zero at all positions dehned by elements in Q. Then we can delete 
those t rows of F, dehned by the elements in Q, and the corresponding columns, as F is 
symmetric, and the corresponding elements of w,w' to obtain and work with F and v,v', 
where v,v' G Lb- Due to the unique mapping from u,u' —)■ —)■ v,v', the commutative 

subgroups from F are in one-to-one correspondence with those from F, but 2* times smaller 
due to the restriction to 0 for t positions, specihed by Q, in u,u'. □ 

So 2* copies of a commuting subgroup of S'-*- comprise a commuting subgroup of S-*-. We 
are particularly interested in the maximum size commuting subgroups of S-*- and S~^, i.e. 
in those subgroups that cannot be contained in larger commuting subgroups. In terms of 
the binary representation, this translates into Ending matrices B where h is maximised, and 
such that BTB^ = 0. Thus our question simplihes to Ending the largest hx n — t B matrices 

such that BTB^ = 0. Remember that e = 


Theorem 2 All maximum commuting subgroups of S-*- and S-*- are of size 2® and 7A 
respectively. Moreover there are Xe -|- 1) such groups in both cases, and each 

element is in n;=;( 2 ^' -|- 1) such groups. 
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Proof. Let the rows, Bj, of B generate the linear subspace, L. Observe that L has a zero 
diagonal. Therefore any Bj G satishes BjTBj = 0, and therefore so does any member 
of L. The maximum size of L is then derived as follows. 

• Choose non-zero Bq G such that BqTBq = 0. There are 2”“^ — 1 such choices. 

• Choose non-zero Bi G F2~^, Bi ^ Bq, such that any member, h, of the linear space 

generated by {Bq, Bi} satishes feTfo^ = 0. There are — 2^ such choices. 

• Choose non-zero B2 G i ?2 not in the linear space generated by {Bq,Bi}, such 

that any member, b, of the linear space generated by [Bq, Bi, B2} satishes bTb^ = 0. 

There are 2”'“* — 2^ such choices. 

«... and so on. 


We continue in this manner until there are no more choices. Thus B always has a maximum 
of e rows, i.e. a maximum commuting subgroup of S'-*- is always size 2®. In total we generate 
11^=0— 2^) linear subspaces. But, for M a maximum rank e x e binary matrix, 
B' = MB generates the same linear subspace as B. We have to remove such repetitions. 
M is chosen from GL(e, 2) and |GL(e,2)| = YYj2o{‘^^ ~ 2"'), so the total number of unique 

linear subspaces, i.e. of maximum commuting subgroups of S-^, is Xe = ^ = 

11^=1 (2“* + !)• For F-*- we simply multiply the maximum commutative subgroup size of S-*- by 
2*, as there are t redundant rows/columns, to get 2”“*^ = 2® x 2h The number of maximum 
commutative subgroups remains at Xe = nk(2^ + i). 

Each element is in 11^=1 (2"' * ^ — 2-') linear subspaces, since we hx the hrst element. 
There are ny=i(2^ ~ 2-^) repetitions (the elements in GL(e, 2) with hxed hrst row), so each 


element is in 


ne-l{2n-t-i_2q 


n-=l(2'=-2q 


= n^=i(2-*' + 1) such groups. 


□ 


It is interesting to note that Xe is also the total number of binary self-dual codes of length 
2(e -|- 1). This leads us to make a passing observation regarding the multiplicative order of 

f. 


Lemma 11 Let u he the multiplicative order ofT, i.e. let r“ = / for some minimum positive 
u. Then u is even and T cannot he factored as T = 1212^ for some 12. 

Proof. If u is odd or if 312 such that T = 1212^ then we show that 35 matrices of size 
xn — t such that BTB^ = 0, where the set of B matrices is in one-to-one correspondence 
with the set of self-dual binary codes of length 2e and dimension e, implying that there are 
Xe-i of them. But this is impossible as, by theorem]^ there are Xe of them. The argument 
is as follows. If u is odd then BTB^ = ct'^ ft= CC^ = 0, where B = Cf'^. If 
T = 1212*^ then BVB^ = CC^ = 0, where C = BVL. In both cases C is taken from the set 
of matrices that generate self-dual binary codes of length n — t and dimension where 
B = Cf in the hrst case and B = 1712“^ in the second case. There are Xe-i such matrices, 
which contradicts theorem [2l □ 
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Example 4 Let G be the graph defining a mixed graph state, p, stabilised by the rows of 

/ X ® z ® z ® z \ 

A= ■ We want to find all maximal commutative subgroups of S-^. The 

\ I®I®Z®X 


adjacency matrix of Gb is T = 


, which has full rank, so F = F. Therefore 


/ 0 1 1 1 \ 

10 10 
110 0 
\ 1 0 0 0 / 

t = 0, e = _ nm _ 2 ^ g, maximal commutative subgroup of S'-*- is 

2 n-e _ 22 _ 4 ^ Moreover, F“ = F^ = J. 

From theorem there are Xe = 15 maximum commuting subgroups of S'-*- and these can 
be generated by 15 matrices, B, satisfying BTB^ = 0, where the B can be chosen, (non- 
uniguely) from 


1000 

0110 

0100 

1011 

1010 

1001 


1000 

0101 

1100 

0010 

1010 

1101 


1000 

0011 

1100 

1001 

0110 

0001 


0100 

1010 

0010 

0001 

0110 

1001 


0100 

0001 

0010 

1101 

1011 

0111 


1100 

1001 


For instance, B = 
maximum commuting subgroup with elements {I ® I 


acting multiplicatively on the rows of generates the 

I, Tiy<ZX®Z®I, =F 


iY ® I ® Z ® X, I ® X ® I ® X}. Similarly, B = 


® I C 

1011 

0111 


acting multiplicatively on the 


rows of A^ generates the maximum commuting subgroup with elements {I ® I ® I ® I, =p 
iX ® I ® Y ® Y, ^ iZ ® X ® X ® Y, ^ iY ® X ® Z ® 1}. Observe that, in this case, 
^iY ® X ® Z ® I occurs in both subgroups. More generally, each group element occurs in 3 
maximal commutative subgroups. 


Example 5 Let G be the graph defining a mixed graph state, p, stabilised by the rows of 


>1 = 



0 1 1 

The adjacency matrix of Gb is T = \ 1 0 1 

1 1 0 


, which has rank 2. 


Therefore t = 1, and e = ^ = 1, and the size of a maximal commutative subgroup of S'-^ is 
2 n-e _ 22 _ 4 ^ Qy removing one dependent row and corresponding column of T (we choose 

0 1 
10 

From theorem^there are Xe = 3 maximum commuting subgroups of and these can be 
generated by extensions of 3 matrices, B, satisfying BTB^ = 0, where the B are 


the last row/column) we can obtain at = 


( 10 ) , ( 01 ) , ( 11 ) 
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Noting that the last row/column ofT is dependent beeause ro + Fi + r 2 = 000, we obtain 
the three maximum commuting subgroups of via the 3 matrices, B', these being 


100 \ 
111 ) ’ 



110 

111 


leading to the following three maximum commuting subgroups of S^, respectively: {J®/® 
I,X^I^Z,TiY^IY0Y,TiZ0Y®X}, {I , Z ®X ,TiY 0IY ^Y,±iX 0 Z 0Y}, 
{I (ZI(ZI,tY (ZX (ZZ, TiY (Z lY (ZY,I (Z Z (Z X}. 

Lemma 12 Let G and H be two mixed graphs on n variables. If they have the same mixed 
rank e, there exists an isomorphism between their commutative subgroups. 


Proof. Assume first n = 2e. Then, the linear group generated by their respective F is 
in both cases F^. The isomorphism is given by sending the basis elements of (F*^) (i.e., the 
rows of F*^) to the corresponding elements in (F'^). Their commntative properties are the 
same, so this isomorphism preserves commntatitivity. 

For n > 2e, the linearly dependent elements can be expressed as F^ = YIk 
similarly (after a possible reordering, that is, a permutation of the nodes) = YliK'^k ■ 
The isomorphism for the independent rows is the same as before. For the dependent rows, we 
send AjAf^ to AjAj^,. This sends the row 00... 0 to itself, and commutativity is preserved 
by linearity. □ 


Example 6 Let Gf, be the line in 4 variables, with adjacency matrix F = 


/ 0 

1 

0 

0 \ 

1 

0 

1 

0 

0 

1 

0 

1 

^0 

0 

1 



and the complete graph in 4 variables. Both of them have e = 2. The group generated by the 

respective adjacency matrices is Fl ITe have the isomorphism f ■. S'^ ■ —)■ Sf-i , , , 

J ^ r j Line Complete graph 

defined by f{Al) = ^{o, 2 , 3 }> fiAf) = fiA^) = ^{ 0 , 1 , 3 }■ particu¬ 

lar, the commutative subgroup on the line generated by Aq,AJ corresponds to the commuta¬ 
tive subgroup on the complete graph generated by A^q 2 ^}^ 


Example 7 let G be the triangle on 3 variables, and H the star graph on 3 variables. For 
both of them, e = 1. A possible isomorphism is given by f : Sq ^ Sjj, defined by /(.Ag) = 
Aq, fiAf) = Af, f{A^) = ‘A|"g 2 }- In particular, the commutative subgroup on G generated 
by Af,A^Q 2 } corresponds to the commutative subgroup on H generated by Af,A 2 - Note 
that if we take the isomorphism given by f : Sq ^ Sjj, defined by /(Alg) = Aq, f{Af) = 
A 2 , /(A,^) = Ai^QiY we obtain the same commutative subgroups. 


5 Children of pure graph state parents 

In this section, given a binary string of length n, j = jo - ■ ■ jn-i, we dehne Sj := Af^, where 
K = {i '■ ji = 1}) where S0...00 = I- Remember that, for i,i' C K, i' > i, we dehne Aj^ to 
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contain the multiplicative factor However the ordering is just as valid and 

differs by a global multiplicative factor of —1 if and anti-commute. This situation 
carries over to the *Sj. In theorem]^ and for such cases where *Sj contains anti-commuting 
factors, the bj are therefore only dehned up to ±1. Let Lm = {n : {A^)v,n+m = Z or Y}, and 
let J = {x e F” I I{x) = 1}, where I{x) = nmio(l + = YlA=o^rn{x). 

Theorem 3 Given a child, p, of a pure graph state parent of a mixed graph, corresponding 
to a symmetric extension of A by e columns and rows, A'^, we can write p as: 


P 


1 

2n 


3&J 


bj e {±l,±*},Vj, 


where all *Sj, j G J, commute pairwise, and if*Sj,*Sj/ are present in the sum, so is *Sj%ji = 
Sj+y, implying that {sj : j G J} is a commutative subgroup of S~^. 


Let G F 2 and G F 2 , 0 < a, (3 < n he weight-one and weight-two binary vectors 
with I’s only in positions a, and a and ft, respectively. 


Lemma 13 Given a child, p, of a pure graph state parent of a mixed graph, described by 
A^, the coefficients, bj, in the sum of theorem^ are as follows: 

• Gase e = 1: 


— If *Sj<y is present in the sum then bj^ = -|-1. 

— If *Sja,f} is present in the sum then we have two cases: 

1. if*Sjc and ^jP anticommute, then 3 a, 6 G {a, (3} such that {A^)a,n = Z and 
(A^)b,n = Y; furthermore, the corresponding term in the sum will be egual to 

= iAfAf, so bjc,p = ±i. 

2. if Aff and Af^ commute, the corresponding term in the sum will be egual to 

— Any other matrices, Aj^, will be products of Aj^. also present in the sum, with Kj 
of size 1 and size 2. Furthermore, their coefficients are given by the multiplication 
of the coefficients of the terms into which it is decomposed. 

• General e: The coefficient of Af^, with K of size 1 or 2 is given by the multiplication 
together of the coefficients, one from each extension column. In this case, there might he 
nondecomposable terms. However, the coefficient of a term will be giving by multiplying 
the coefficients obtained by taking each column separately. 

• Furthermore, for any e, if we change the sign of a term Af, 0 < k < n, then we add 
the Boolean linear term Xk to the quadratic Boolean representation of the pure parent 
graph state, or if we change the sign of Af^, for K a set of size t > 1 present in the 
sum, we add any of the Boolean linear terms Xk for each k s.t. {A^)n+j,k = Y or 
{A'^)n+j,k = Z. 
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Example 8 Let the direeted triangle he defined by the stabilizer basis AQ = X®Z®fiAi = 
I ® X ® Z, A 2 = Z ® I ® X. Then the basis of S~^ is given by reversing the arrows: 
Al = X ® I ® Z, Af = Z ® X ® fi Al = I ® Z ® X, so = {sqoo = / ® ® /, sioo = 
X®IZ>Z, soio = Z®X®fi siio = -iYZ>X®Z, sqoi = I®ZZ>X, sioi = iX®Z®Y, son = 
-iZ ®Y ®X, Sill = -iY O F 0 F}. 

One of the two parents is formed by adding the column (X Z F)^ to A, giving the parent 
j2{xoX2+XiX2+XiX3+X2X3+X2)+X2 _ j 2 {xf,X 2 +XlX 2 +X 2 + (X{x) + l)xz)+X 2 /.g = |1 2} and L[x) = 

xi + X 2 + I, so J =< 100,011 >= {000,100,011, 111}. IFe can re-interpret I{x) and J as 

parity and generator matrices, H and G, respectively, where H = { Oil ) and G = ^ 

and where G generates the binary linear code with codewords in set J. By tracing over the 
environmental guhit, X 3 , we get 


p — — + + Z) X Y ZY Z F) — — (sgoo F '^loo F isgn F 


Note that XZiZZ and ZzYcommute, and that {XZIZZ) ■ {ZZYZX) =YzYZY. 
We can obtain another parent for the same child by doing a local complementation (LC) on 
vertex 3 (i.e. LC 3 ) on the graph described by the parent. This is an LC in the environment. 

If one describes the Tj^-linear offsets of the parent by black vertices in the graph, then LCs 
swaps the vertex neighbours of vertex 3 from black <H- white. 

The other of the two parents is formed by adding the column (X F ZY to A, giving 
the parent jfiOox 2 +xi+xixz+x 2 xz)+xi _ j 2 {xox 2 +xz+{T(x)+i)xz)+xi ^ Once again L 3 = {1,2}, T{x) = 

Xi+ X 2 + 1, and .J = {000,100,011, 111}. By tracing over the environmental gubit, X 3 , we 
get 

p=\{I Zl Zl Y X Zl Z Z - Z ZY ZX - Y ZY ZY) = \ (sqoo F sioo - isoii - ism) . 

o o 


H and G are unchanged. Observe that the second parent is obtained from the first by swapping 
the positions of Y and Z in rows 1 and 2 of the extra column (column 3). In terms of 
the Boolean function representation of the parents, this translates to adding the Boolean 
guadratic term X 1 X 2 + X 1 + X 2 , and removing X 2 and adding Xi Z^-linear terms. Observe that 
the coefficients of the Pauli basis terms Z ZY Z X and Y Z Y Z Y are multiplied by —1. 
This is because both terms are generated from rows 1 and 2 of A^ which are the rows where 
Y and Z are swapped in A^. 

For each parent we can also consider how the addition of binary linear terms (in the lab) 
affects the resultant child density matrix. We don’t currently consider the addition of binary 
linear terms in the environment (i.e. X 3 for this example). For instance, for the addition of 
column (X Z F)^ j 2 (xox 2 +xiX 2 +xiX 3 +x 2 X 3 +xo+x 2 )+x 2 ^ addition of Xq, simply flips 

the signs of X ZIZ Z and Y ZY ZY as both terms have row 0 of AA as a factor. However 
addition of xi or X 2 has the same effect as swapping Y and Z in rows 1 and 2, so swaps 
between the two parents. Thus the addition ofX{x) = X{x) + 1 fixes the child density matrix. 
So we have the following maps for a G {0,1}.' 
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child = 8/9 


parent 


1 '^' _ 

2{xoX2 + X 1 X 2 + X 2 + I{x){x 3 + a)) + X 2 , 

2(xoX 2 +Xi+ X 2 + X(x)(x 3 + a)) + Xi, 

2(xoX2 + X 1 X 2 + Xo + X 2 + X(x)(x 3 + o)) + X 2 , 
2(xoX2 + Xo + Xi + X 2 +X(x)(x 3 + a)) + Xi, 

2(xoX 2 + Xi + X(x)(x 3 + a)) + xi, 

2(xoX2 + X 1 X 2 + Xi + X 2 + X(x)(x 3 + a)) + X 2 , 
2(xoX2 + Xo + Xi + X(x)(x 3 + a)) + xi, 

2(xoX2 + X 1 X 2 + Xo + Xi + X 2 + X(x)(x 3 + a)) + X 2 . 


I ® I ® I + X ® I ® Z + Z ®Y ® X + Y ®Y ®Y 


I ® I ® I - X ® I ® Z + Z ®Y ® X - Y ®Y ®Y 


I ® I ® I + X ® I ® Z - Z ®Y ® X - Y ®Y ®Y 


I ® I ® I - X ® I ® Z - Z ®Y ® X + Y ®Y ®Y 


Other children given by the extensions {Z X Y)'^ and {Z Y X)^ are respectively 


Pi = i(^ 
P^ = W 


I Z) I + aZ Z) X Z) I + hX Z) Z Z)Y + cY Z)Y Z)Y) = \ (sooo + asow — biswi + cisni) 
I Z d Y aY Z X Z Z Y bl Z Z Z X Y cY ZY ZY') = 000 “ 1 “ -s 110 “ 1 “ ^"Sooi Y ci s 

with condition c = ab, a, 6 G {1, — 1} in both cases. 


( 


Example 9 Let a mixed 6 -clique graph be defined by the stabilizer basis A = 


V 


X 

I 

I 

I 

I 

I 


z 

X 

I 

I 

I 

I 


z 

z 

X 

I 

I 

I 


z 

z 

z 

X 

I 

I 


z 

z 

z 

z 

X 

I 


z 

z 

z 

z 

z 

X. 


\ 


Then the basis of is obtained from . The parents are obtained by adding e = 3 columns 
to A, with subsequent addition of e = 3 rows. For instance, one parent is given by 


A^ = 


/ X 

z 

Z 

Z 

z 

z 

X 

I 

X \ 


/ X 

I 

I 

I 

I 

I 

I 

I 

I 

\ 

/ 

Y 

Z 

z 

z 

z 

Y 

I 

X 


I 

Y 

I 

I 

I 

I 

Z 

I 

I 


I 

I 

X 

z 

z 

z 

z 

I 

X 


I 

I 

X 

I 

I 

I 

Z 

I 

I 


I 

I 

I 

X 

z 

z 

I 

X 

Y 


I 

I 

I 

Y 

Z 

Z 

I 

I 

z 


I 

I 

I 

I 

X 

z 

I 

Y 

Y 

= 

I 

I 

I 

Z 

X 

z 

I 

Z 

z 


I 

I 

I 

I 

I 

X 

I 

z 

Y 


I 

I 

I 

Z 

Z 

Y 

I 

z 

z 


I 

Z 

Z 

I 

I 

I 

X 

I 

I 


I 

Z 

Z 

I 

I 

I 

X 

I 

I 


I 

I 

I 

I 

z 

z 

I 

X 

I 


I 

I 

I 

I 

Z 

Z 

I 

X 

I 


\ I 

I 

I 

z 

z 

z 

I 

I 

X i 


\ I 

I 

I 

Z 

Z 

Z 

I 

I 

X. 

/ 


which represents fiYi^&+^2X6+x3X4+x3x^+x3xs+x4X5+x4X'j+x4xs+x5XT+x5xs+xi+xs+x4+x5)+xi+x3+x5 Here 
L3 = { 1 , 2 }, L4 = ( 4 , 5 ), andL^ = { 3 , 4 , 5 }, andX{x) = (xi+X2+l)(a;4+a;5+l)(x3+X4+X5 + 

1), so J =< 100000, 011000, 000011 >= {000000,100000, 011000,111000,000011,100011, 011011,111011}. 
We can re-interpret I{x) and J as parity and generator matrices, H and G, respectively, 

011000 \ / 100000 

where H = | 000011 j and G = j 011000 | , and where G generates the binary linear 
000111 / \ 000011 
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code with codewords in set J. By tracing over the environmental qubits, Xq,x-j,Xs, we get 

p= {I ® I ® I ® I ® I ® I + X ® Z ® Z ® Z ® Z ® Z - I ® X ®Y ® I ® I ® I 

-X ®Y ® X ® Z ® Z ® Z - I ® I ® I ® I ® X ®Y - X ® Z ® Z ® Z ®Y ® X 
+I ® X ®Y ® I ® X ®Y + X ®Y ® X ® Z ®Y ® X). 

Proof. (Theorem [^: The qubits present in each Lm will be the ones connected with 
the environmental qubit n + m, implying that the difference between the measurement in 
Xn+m = 0 and Xn+m = 1 will be the Boolean function Yhk&Lm support of p will be 

equal to the binary vectors such that these Boolean functions are 0 for all m. By lemma 
p can be expressed as the sum of some of the matrices on S'-*-; by lemma all the matrices 
in S-*- have nonintersecting support: p will be equal to the sum of the matrices in whose 
support intersects the support of p; therefore, it is sufficient to find the nonzero entries of 
the first row of the matrix. It is easy to check that any matrix in has only one nonzero 
entry on the first row. In other words, p will be equal to the sum of the matrices whose 
support intersects {(0, k), k G K}, where K = {X]m=o ^ J}^ and J 

e—1 

is the indicator of the Boolean function n(E + 1). To prove that these matrices are 

m=0 k£Lm 

the Sj, where j G J, we need the following lemma: 

Lemma 14 Let sj G S-^. Then, {sj)o,t = l ^ k = for j = joji . ..jn-i- 

Proof. Let *Sj = <S)rx ^ ^Ry ^ ^Rz ^ ^Ri M he a. t x t matrix. Then, 

^ f M 0 \ ^ f 0 M \ 

0 M ) 0 j 


Note that X has the same support as Y, and Z has the same support as /, so we only have 
to consider this two cases. This implies that support(M 0 J) =support(M), while (0,/c + 
2*“^) Gsupport(M (g) X) (0, fc) Gsupport(M). This implies that {0,k) Gsupport(sj) 

k = J2 Z&RxURy 2™"^ = E™=0 for J = JoJl ■ ■ • Jn-l- 

□ 

To complete the proof of theorem]^ there remains to prove that all Sj, j ^ J commute 
pairwise, and if Sj, are present in the sum, so is 

First, we shall prove that if s^, are present in the sum, so is this follows from the 

indicator J being a linear space, since corresponds to the sum of the Boolean vectors. 

Now we shall see that it is commutative: Let j,k E J: We can write Sj = and 

*Sk = riftes where Af, Aj are in the basis of S'-*-. Then, the corresponding Sla G S, a G v4, 
have an even number of Y and Z in each extension column, and similarly for B. Therefore, 
both *Sj and have either X or / in all extension columns, and therefore they commute. 

□ 


Proof. [Lemma 13 Since all the parent Boolean functions have no constant terms and no 
linear terms involving the environmental qubits, the first entry in a truth table for each mea¬ 
surement Ifjm) will always be -|-1. The first row of the density matrix for each measurement 
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will therefore be equal to Also, note that the hnal p will be nonzero only 

where the entries for all \'iprn) are equal, and will then be equal to any of them; it is therefore 
enough to look at ('0o...o|- 

Case e = 1; 

• Size 1: By the proof of theorem any term, 0 < j < n — 1 is such that 

{A^)j^n = X OT I. Since row j of the stabilizer of the parent graph state, is equal to 

(extension entries), there is no linear term xj in the parent graph, so that the entry 
0... 1... 0, where the 1 occurs in position u, of |' 0 o...o )5 and therefore of ('0o...o| of the 
parent is As the hrst row of Aj has in the same position, the coefficient 

of Sj is +A. 

• Size 2: Suppose that the matrix Ajjj^y 0 < j, A; < n — 1, is present in the sum. Then; 

— if Aj and A'^ anti-commute, the term in the sum will be 0 < j, fc < n— 1 , 

since otherwise ±^AjA^ 0 < j, fc < n — 1 would not be Hermitian. Furthermore, 
neither Aj or Aj are present in the sum, because by theorem this would imply 
that both would be present and would therefore be commuting. This implies that 
3 a, 6 G {j, k} such that {A^)a,n = Y and {A^)b,n = Z. The hrst row of ±^AjA^ 
has its only nonzero entry where Xj = x*, = 1, = 0 Vu 7 ^ j,k- Therefore, 

the entry in will be given by the presence or absence of XjXk^Xj^Xk, in the 
ANF of the parent graph. Since Aj and A'^ anti-commute, the restriction of 
Aj,Ak (that also anti-commute) to the pair j,k will give XaZj^ and laXp for 
a, f3 E {j,k}. This implies that the restriction of A'^ and A'p are , respectively, 
Xalp and ZaXp, so the restriction of iA'^A'^ is Y^Xp. Note that we only need to 
look at the restriction, since any further Kronecker product will be by X, Z or I, 
and these will not change the hrst non-zero entry. 

* If the extension column gives (A^)a,n = Z and (^'^)/ 3 ,n = Y, we get the term 
2 {xaXi 3 + xp) -|- Xjs in the parent graph, so |V’o) has entry i for Xj = Xk = 
l,Xu = 0 yu ^ j, /c, which implies that (- 00 1 has entry —i, same as the same 
entry in the hrst row of the matrix given by the order of multiplication A^A^. 

* If the extension column gives (A^)a,n = Y and (^®)/ 3 ,n = Z, we get the term 
2xa + Xa in the parent graph, so l-^o) has entry —i (so {'ipo\ has entry i) for 
Xj = Xk = l,Xu = 0 Wu ^ j, fc,same as the same entry in the hrst row of the 
matrix given by the order of multiplication A^A'^. 

— Suppose Aj and Ak commute. Then, the term < j, fc < n — 1 is 

present in the sum, the matrices are commuting since AjA^ is Hermitian. By 
inspection, on the pair Aj,Ak E S, X^Zf^ and Z^Xfi give entry -1 regardless of 
extension (note that they have the commuting entry in the extension column), 
same as AjAk. As for Xaljs and laXjS give entry -|-1 regardless of extension, 
same as Ae^A^- 

— By the proof of theorem]^ A"^, K C {0,..., n — 1}, is present in the sum ih the 
corresponding stabilizer, ^^,has either X or / in the extension column. Therefore, 
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there exists a (not necessarely unique) decomposition in size 1 and size 2 terms 
that have either X or / in the extension column (since YZ = iX), and, as the 
entry xj = 1 Vj G A, Xj = 0 Vj ^ A will depend on the sum of the terms for 
the size 1 and size 2 cases, as the Boolean function has degree at most 2, this 
entry will depend on the entry of the size 1 and size 2 cases, and therefore the 
coefficient will be the multiplication of the size 1 and size 2 cases (for any given 
decomposition). 

General e: Since the case e = 1 was independent of the actual stabilizers (it only depends on 
the extension), each new column will modify the Boolean expression accordingly. Therefore, 
the coefficient of any term is given by the multiplication of the coefficient resultant of each 
extension column. 

Change of sign: Let us first assume that j has weight 1, with support in k. Note that 
each Aj has support in ( 0 , 2 ^). In terms of (t/’o...o|; this means that changing the sign of Aj 
changes the sign in the jth element in |t/’o...o) (and therefore in (t/’o...o|)- This is equivalent to 
adding a linear term Xk as long as we change the sign all elements in|'^o...o) (and therefore 
in (t/’o...o|) with Xk = 1- Suppose now that j has weight > 1. If we change the sign in j, we 
have to change also the signs of an odd number of its decomposing terms. Consider therefore 
the smallest terms present in the group. If the sign changed is of weight 1, see above. If 
the weight is t > 1 , with support in ki,... ,kt, then any of the Aki are not elements in the 
subgroup. Adding any linear term Xk^ will give the desired matrix, since it is zero in the 
places were it might differ. 

□ 


Corollary 1 The commutative subgroup corresponding to a child is maximal. 

Proof. Each column added gives a affine linear Boolean function /j, which gives a con¬ 
strain to the indicator J of f = Yi ft- Each constraint, /j, is nontrivial (that is, not a 
constant), because if this were the case, then the column would only have X and I, and 
would therefore be superfluous, yielding a contradiction. The fi are all independent, other¬ 
wise we get redundant columns, yielding a contradiction with e being minimal. Furthermore, 
any new independent linear constraint reduces by half the size of the indicator, wich means 
that this size is equal to 2 ”“®, so it is a maximal commutative subgroup. □ 

NB: Note that allowing superflous constraints fi will give commutative subgroups that 
are in general not maximal. In this way, we could also extend graph states in a natu¬ 
ral way: the density matrix for the graph state is given by the sum of all the elements 
of the stabilizer (since it is self-dual), and of course any consistent sign changes that give 
linear terms (eigenvalue -1). We can however dehne more density matrices that are sta¬ 
bilized by the stabilizer of the graph state, by allowing also smaller commutative sub¬ 
groups. For instance, a density matrix associated to the undirected line from 0 to 1 could be 
then p = ao{I<Zl±X<ZZ±Z<ZX±Y<ZiY) + ai{I<Zl±X<ZZ) + a 2 {I<ZiIdxZ<ZiX) + 
a^{I I ±Y ® F) -|- 04 / 0 /, where ^ a* = 1. Note that in the hrst term, the signs have to 
be consistent, changing for instance X Z) Z to —X Z) Z forces the change Y Z)Y to —Y ® Y. 
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6 Open problem: the weighted sum of a maximal com¬ 
mutative subgroup of S~^ is a density matrix 


Conjecture 1 Any maximal commutative subgroup corresponds to a child of a pure parent 
graph, so its weighted sum (with appropiate coefficients) is a density matrix. 


Froo/;[incomplete] Let be a maximal commutative subgroup of for a mixed graph 
G. Then, the corresponding elements in S commute as well (because direction of arrows is 
not important), and form a commutative subgroup, M. If there is an extension of S (not 
unique) to the stabilizer of to a pure graph state, the extension has e columns. Since the 
elements of M commute, for e = 1, we can assign X for all odd products of elements of the 
basis S, and / to even products. Since it is a subgroup, this is consistent. In general, we 
can at least always assign either X or J in the corresponding columns of the basis, which 
respects their commutativity. Not all combinations will be valid (for instance, assigning 
to a node that does not commute with everybody is not allowed) but there exist at least 
one such combination, by lemma In general, we still have to prove that this combination 
is always possible.If it exists, any such combination would have as density matrix the one 
generated by M^. 

A possible strategy for hnding these extension columns is to write the indicator that will 
indicate the position of the Y, Z positions. For any nondecomposable Jff), we need to impose 
the condition that in the indicator function, all the Xj, j G V, are equal. This we can achieve 
by taking the products n(^i + + !)• This strategy works for non-intersecting elements. 


Lemma 15 Regardless of whether we get intersecting elements or not, we always obtain e 
terms in the final product of the indicator. 

Proof: Let M be a maximal commutative subgroup. By theorem \M\ = 2"^“®. Let C 
be the corresponding binary linear code. A binary linear code with k = n — e independent 
codewords has a parity-check matrix with e rows: the product of any set of these e parity- 
check conditions will be our indicator function. 

Note that since e < | by dehnition, we have 2e < n, which means that e < n — e. 


Example 10 Consider the mixed graph given by 


/ X®Z®I®I®Z 

— -Ao \ 

I®X®Z®I®Z 

= Ali 

I®I®X®Z®Z 

= C\.2 

I®I®I®X®Z 

= ^3 

\ I ®I ®I ®I ®x 

= -Ai / 


Then, the dual is 


/ X ® I ® I ® I ® I 


Z ®X ® I ® I ®I 

= Jif 

I®Z®X®I®I 


I®I®Z®X®I 


\z®z®z®z®x 

= AJ J 
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The rank of the corresponding adjacency matrix is so e = 2 . Commuting with Aq is 
the subgroup generated by A^, AJf, A^, Af^i 4 ^ = Ae^Arf. But not all of them commute so we 
choose another element in this subgroup. Commuting (e.g.) with both Aq and Aj is the 
subgroup generated by , ^3 , 2 4 }- ® maximal commutative subgroup. We shall 

find an extension of the stabilizers such that we obtain a parent pure graph state. 

For any nondecomposable Affj, we need to impose the condition that in the indicator 
function, all the Xj, j G D, are egual. This we can achieve by taking the products n(^i + 
kk + l). Here we have the composite Af^^ 24 }; which gives the condition (a;i+a ;4 + l)(a; 2 +a ;4 + l) 
(or any of the eguivalent conditions). We assign then in the first column of the extension the 
entries Z, Y to positions 1,4, while the remaining entries of the first column will be X or I. 
Since rows 1 and 4 anti-commute, we assign Z to position 1 and Y to position 4 (note that 
the other noncommutative choice will give the same commutative subgroup, with a change 
of sign in ^|'i 2 4 }/ Both Aq and A 2 anti-commute with A(^ and Aj, so we assign X to 
positions 0 and 2. Finally, we assign I to position 3. Now for the second column: Here we 
assign Z or Y to positions 2 and 4, while we assign X or I to the remaining nodes. Since A '2 
and A'l commute, we give them the same entry, say Z. Now, A'q is already commutative 
with all the other rows, so we assign I to position 0, and I to position 1 for the same reason. 
We assign X to position 3, since A\!q anti-commutes with both Ai 2 (ind Aj^ . This is now a 
fully commuting set: 


/ X®Z®I®I®Z 

X I \ 

I®X®Z®I®Z 

Z I 

i®i®x®z®z 

X z 

I®I®I®X®Z 

I X 

\ I ®i ®i ®x 

Y Z ) 


Example 11 Let G be the arrrowed 6 -cligue given by the following stabilizer basis: 

/ X z z z z z \ 

I X z z z z 

I I X z z z 

I I I X z z 

I I I I X z 

\ I I I I I X ) 

Then, e = 3 . ^4 maximal commutative subgroup is generated by .Ag , 2}) *^{3 4}- 

the indicator {xi + X2 + l)(a^3 + X4 + l)(a;5 + 1 ). According to this, we assign in the first 
extension column the entries Z and Y to positions 1 and 2, respectively since their respective 
rows anti-commute. The remaining entries on this column have to be either X or I. We 
assign X to positions 0,3,4,5, because their respective rows anti-commute with both. Now 
for the second column: We assign Z and Y to positions 3 and 4, respectively since their 
respective rows anti-commute. The remaining entries on this column have to be either X or 
L We assign X to positions 0 and 5, because their respective rows (still) anti-commute with 
both, I to positions 1 and 2, which respective rows now commute with all rows. Finally, in 
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the last column we assign Z to position 5. The remaining entries on this column have to he 
either X or I . We assign X to position 0, and I to all the other positions, since their rows 
all commute pairwise now. This is now a fully commuting set. 


/ 

X 

z 

z 

z 

z 

z 

X 

X 

X 

\ 


I 

X 

z 

z 

z 

z 

z 

I 

I 



I 

I 

X 

z 

z 

z 

Y 

I 

I 



I 

I 

I 

X 

z 

z 

X 

z 

I 



I 

I 

I 

I 

X 

z 

X 

Y 

I 


V 

I 

I 

I 

I 

I 

X 

X 

X 

z 

) 


7 Appendix A 


We re-iterate that, as with S, if some or all the rows of A pairwise anti-commute, then 
some or all of the members of S'-*- are dehned up to a global constant of ±1. This does not 
affect the commutation relations between members of S and S'-*-. But which member^of S'-*- 
are dehned up to ±1 and how many of them? Consider the n-vertex mixed graph G with 
adjacency matrix , derived from A^. Then Gb is also the undirected graph associated 
with G (as well as with G), as Gb has adjacency matrix T = A -|- A^. Then the members 
of S'-*- associated with a -1-1 coefficient are obtained as the subset of row products of A'^ 
indexed by sets representing all possible independent sets in Gb- But how to compute these 
sets? For a set w, the power set of w, V{w) is the set of all subsets of w including the empty 
set and w itself. Dehne C to be a family of subsets of specihcally let V represent the 
maximum independent sets in Gb meaning that, for v C Zn = {0, 1,..., n — 1}, v G V iff 
V is an independent set in Gb and there does not exist a v' which is an independent set in 
Gb such that v <Zv', v ^ v'. Then Gb specihes V precisely and V specihes Gb precisely. For 

/ 0 0 0 1 1 
0 0 0 0 1 


instance, for n = 6 let Gb be dehned by the adjacency matrix F = 


V 


1 \ 
0 
0 
0 
1 
0 / 


Then V = {{0,1,2}, {2, 3,4}, {1, 2,3, 5}}. Observe that V also lists the maximum cliques in 
the graph complement, Gb of Gb- From V we can compute those row subsets of A^ whose 
product is independent of product order, i.e. the row subsets that are fully commutative. 
We call this set of subsets S(V) and a member of S(V) represents a unique member of S'-*- 
with a -|-1 coefficient. For our example there are |T(0)| = 24 such row subsets, namely 

g(V) = {0,{O},{1},{O,1},{2},{O,2},{1,2},{O,1,2}, 

{3}, {1,3}, {2, 3}, {1, 2, 3}, {4}, {2,4}, {3,4}, {2, 3,4}, 

{5}, {1,5}, {2, 5}, {1,2, 5}, {3,5}, {1,3, 5}, {2, 3, 5}, {1,2, 3, 5}}. 


Then S(y) = 'P(Z„) \ S(y) is the family of subsets of Z„ that represent non-commuting 
row subsets of AA, i.e. a member of £{y) represents a unique member of S'-*- with a ±1 
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coefficient. The number of them is |£^(V^)| = 2"' — |T(i/)|, this being the number of members 
of S~^ with a ±1 coefficient. For our example |T(1/)| = 2® — 24 = 40. We can express S in 
terms of power sets as follows, 

£{V)= \JV{w). 

wGV 

We now provide a recursive algorithm to compute T(14), given (and therefore V). Let 
w GZn,V = \Jw,w'&v,w^w' wnw' and, for w eV, let 14 = ^ example, 

y = {{0,1,2}, {2,3,4}, {l,2,3,5}},l> = {{0,l,2}n{2,3,4},{0,l,2}n{l,2,3,5},{2,3,4}n 
{1,2, 3, 5}} = {{2}, {1,2}, {2, 3}} and = {{0,1, 2} n {2, 3,4}, {0,1, 2} n {1, 2, 3, 5}} = 

{{2}, {1,2}}. Then, recursively, 

£{V)=£{V)U \JV{w)\£{V^). 

wGV 

So, for our example with V = {{0,1, 2}, {2, 3,4}, {1, 2, 3, 5}} we obtain 

= T({{2}, {1, 2}, {2, 3}}) U Vm 1, 2}) \ ^({{2}, {1, 2}}) 

U V{{2, 3,4}) \ T({{2}, {2,3}}) U V{{1, 2, 3, 5}) \ ^({{1, 2}, {2, 3}}) 

= T({{1,2}, {2,3}}) U V{{0, 1,2}) \ V{{1, 2}) 

U V{{2, 3,4}) \ V{{2, 3}) U V{{1, 2, 3,5}) \ T({{1,2}, {2, 3}}) 

= {{0},{0,l},{0,2},{0,l,2}}U{{4},{2,4},{3,4},{2,3,4}}UiP({l,2,3,5}), 

from which we obtain the family of 24 sets written above. 


8 Appendix B 


8.1 e = 1 qubit extensions 

We consider here the simplest case where e = |rank(r) = 1. 


Lemma 16 Let j be a vertex of the mixed graph described by A, for which e = 1, such that 
vertex j is either disconnected or connected to other vertices only via undirected edges. Then 

AV = 


Proof. From the conditions in the lemma, row j of A commutes with all other rows in A. 
As e = 1 then there are at least two other rows of A, rows h and k, such that 4 
Al nj'^kn ^ ■^}- As we require row j of A^ to commute with rows h and k of A^ and, 

as row j of A commutes with rows h and k of A, then the only choice for A^j n is /. □ 


A consequence of Lemma IT is that, given a mixed graph described by A, we may 
immediately, and wlog, set all Af^ = / for unconnected and undirected vertices j and focus 
on solving the extension problem for the residual directed graph comprising vertices involved 
in one or more directed edges. This is implicitly the case for Theorem below. 
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Theorem 4 Let A be an n x n Pauli stabilizer matrix such that, from e = 1. Then 
A can always be extended to a fully-commuting stabilizer matrix, A^ of the form A^^ = 
, Such a matrix stabilizes a graph state of n + e = n + l gubits, and has 
elements from {X,Y} on the diagonal and elements from {I,Z} off it. 


Proof. (of Theorem 

Lemma 17 For A'^ = [A^j G {/, X, Z,Y},0 < i, j < n), we can wlog fix A^^ ^ = X. 

Proof. (of LemmafTTl) If .4.® „ = X then we are done. If 4,® „ G {Z, Y} then, by Lemma 
we can conjugate .4^ to X and, by Lemmathis leaves the density matrix stabilized 
Dy A unchanged, and preserves the pairwise commutation properties. If A^ n = then 
we multiply row n by any other row, j, where Af^ ^ {Zi,Y,Z}, i.e. „ ^ 

At least one such row, j, must exist for e > 1. This operation preserves the codespace of 
A^, and results in A] 
discussed previously. 

Although Lemma 17 restricts 4.® „ = X we emphasise that there is nothing special 
about X here, we could’ve chosen to conjugate to Y or Z. But fixing to X simplifies 
the problem for us without costing anything in terms of generality, as, by Lemma the 
associated density matrix is unchanged. 


n,n 


G {X,Y,Z}. Then we continue with conjugation if necessary, as 

□ 


Lemma 18 If the rows i and j of A are pairwise non-commuting then 




e 

■j,n' 


Proof. (of Lemma 18) If rows i and j of A, i j, do not pairwise commute then, 
if 4®^ = Af^, or if 4f ,j and/or Af^ equals I, then rows i and j of 4® do not pairwise 
commute either. But we require 4® to be fully pairwise commuting. So 4f„ 7 ^ 4®„ and 
AU,Al„e{X,Z,Y}. ' ' □ 


Lemma 19 If the rows i and j of A are pairwise commuting then 


= 4 ® 


or at least one of 4/„ or 4® „ eguals I. 




Proof. (of Lemma 19) Rows i and j of 4® must pairwise commute and, if rows i and j 


of 4 pairwise commute, then matrices 4/„ and Af^ must also pairwise commute, hence 
the conditions of the lemma. □ 


Lemma 20 Wlog, the elements A^j, 0 < j < n can be restricted to elements from {/, Z}. 


26 





Proof. (of Lemma The last row, .4,® = < j < of A^, must pairwise 

commute with each oiuie n previous rows. We have, wlog, already hxed A^^ = X. 
We choose elements of A^ j, 0 < j < n from {/,X, Z, F}. For each A^j G {X,Y} we 
successively replace the nth row of A^, i.e. A^, with the product of rows j and n, i.e. 

Ae V 


which ensures that, now, 4.® • G {/, Z}, and continue in this fashion until 


At , G {I,Z}, 0 < j < n 


This process preserves the pairwise commuting property for 

□ 

The row product operations of lemma 20 may change 4,® „ = X to At „ G {Z, F} in which 


the rows of A^ 


case, using Lemma IT, a further conjugation operation on the nth column of A^ is required 
to change Atn back to X. Note that, by Lemmasuch an operation leaves unchanged the 
density matrix, p, described by Tl. 

At this point A^ has been transformed, by row operations and conjugation of the last 
column, to a matrix with X or F on the diagonal, and with {/, Z} off the diagonal, apart 
from the last column where A^n ^ {-^5 At,Y, Z}, 0 < j < n, and A^n = X- If A^n ^ ^} 

then we replace the jth row thus A^j —>■ A^A^j so that, now, A^n ^ U, Z}. 

At this point A^ has been transformed, by row operations and conjugation of the last 
column, to a matrix with X or F on the diagonal, and with {/, Z} off the diagonal. 

Lemma 21 A^ is of the form A^ = so must be a symmetric matrix. 


Proof. (of Lemma 21) All rows of A^ pairwise commute. Given that 4.® has X or F 
on the diagonal and {/, Z} off the diagonal, then this is only possible if A'^ is symmetric. 
□ 

Combining Lemmas IT, 1^, M, and ‘A. one arrives at the proof of Theorem]^ □ 


8.2 Conditions for G to have mixed rank e = 1 

Remember that G is the mixed graph dehned by adjacency matrix A, and Gb the undirected 
graph with adjacency matrix T = A + A^. 

Lemma 22 G has mixed rank e = 1 ^ Gb is a complete tripartite grapf^ with possible 
isolated vertices, and where one of the partitions can be empty (so that Gb is actually a 
complete bipartite graph). 


Proof. If (i, j) is an edge in Gb, then rows i and j of the stabilizer basis are anti-commuting 
so, by lemma 
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the corresponding elements in the extension column of A, (A.f n,Tl®„) can 
only take the pairs (X, Z), (X, F), (Z, X), (Z, F), (F, X), and (F, Z). In other words the edge 
{i,i) in Gb forces vertices i and j to be associated with different operators in the extension 
column At - one can think of X, Z, F as three colours. Therefore e = 1 is only possible if Gb 
is three-colourable, which implies that the graph is empty (1-coloured), bipartite (2-coloured) 
or tripartite (3-coloured). If G;, = 0 then rank{T) = 0, so we discard this option. Isolated 


^Thanks to Jon Eivind Vatne for useful discussions about graph theory. 
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vertices correspond to the rows that commute with all other rows, and can be extended by 
the column entry I, as seen in lemma 16 Let us now consider the connected part of Gb- 


We are now going to prove that, if G has mixed rank e = 1, then the connected part of 
Gb is complete tripartite or complete bipartite, with possible isolated vertices; since it has to 
be 2- or 3-coloured, we have to prove that all the vertices on each partition are connected to 
all the vertices in the other partition(s). Consider any pair of vertices (m, v) such that u and 
V are in different partitions. The partitions are assigned two different colours/entries for the 
extension column, in the set {Z,X,Y}. Therefore the rows u and v anti-commute, since 
the entries in the extension column anti-commute. Therefore there must be an edge between 
u and V in Gb, which proves that the graph is complete bipartite or complete tripartite. 

We will now prove that if Gb is a complete bipartite or tripartite graph with possible 
isolated vertices, then G has mixed rank e = 1: 

In general, the extension of G depends only on Gb- We have that {u, v) G Gb iff row u 
anti-commutes with row v. This is true because all unitaries Uj that form the stabilizer basis 
for a mixed graph state have exactly one X in position j, and Z and I in the other positions. 
If (m, v) G Gb, we have in rows u and v the pair X^Z^, luXy {u —)■ v) or X^Iv, ZuXy (m r>), 
and since all entries other than u and v commute (the set {Z, /} is commutative), then row 
u anti-commutes with row v. On the other hand, if row u anti-commutes with row v, then 
one of the pairs X^Zy, lyXy or Xyh, ZyXy must be present in rows u and v, because as 
discussed all entries other than u and v commute, and these pairs correspond to m —)■ n and 
u ^ V, so {u,v) G Gb- This means that the choice of the extension column/row will depend 
solely on Gb, since the only factor is the commutativity/anti-commutativity of the rows. 

The extension of G by a node depends on whether we can extend the matrix by one 
column. We assign to each of the partitions the colours Z, X, Y- Since the graph is complete 
bipartite or tripartite, then any two nodes in the same partition get assigned the same 
colour/entry, and any two nodes in a different partition get assigned a different colour/entry. 
In particular, this means that there are no arrowed edges between any two nodes in the same 
partition, which means that their rows commute. Since they have the same extension entry, 
their extended rows also commute. This also means that there is an arrowed edge between 
any two nodes in different partitions, which means that their respective rows anti-commute. 
Since the extension entry is different for each of these two nodes, and is in the set {Z,X, Y}, 
this means that the extension rows now commute. This proves that one extension column is 
enough. 

□ 


8.3 On parents and children for e = 1 

In general, as we shall now show, for e = 1 we obtain (up to local equivalence) a maximum 
of three mixed states, obtained from a maximum of three parents (up to local equivalence), 
iV'e)) as stabilized by A^- 

^Note that I is not present because since it does not change the commutativity it can only be assigned to 
the rows that connect with any other rows (which correspond to isolated nodes on Gt), and we are considering 
here the connected part of Gb- 
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Lemma 23 A maximum of 6 distinct children is possible when e = 1. 


Proof. Isolated nodes in Gf, have the extension I, so we consider only non-isolated nodes. 
If G/, is not connected then there are two or more non-empty subgraphs that are not connected 
to each other; however, then rank(r) > 2, which contradicts e = 1. This is because G/, is a 
nondirected graph for which any non-empty subgraph has rank at least 2, and the rank of 
Gb is the sum of the ranks of the unconnected subgraphs. 


Lemma 22 states that G is at most 3-coIourable if e = 1, as Gj, is complete tripartite (or 


complete bipartite). Gb determines univocally whether two nodes have the same colour or 
not. There are 3 choices for the hrst colour, 2 for the second, and 1 for the remaining colour. 
In total that makes in principle 6 possible extension columns, yielding 6 different parents. 
Note that we only consider choices for A^j, for j = 0,... ,n — 1. The reason for this is that 
that, as shown previously, we can £x 


□ 


Lemma 24 When e = 1 we obtain at most 3 local unitary inequivalent children. 


Denote the resultant matrix after applying permutation [Z Y) on A, 

i/e 


e 


J = 


Proof. 

0 ,..., n — 1, as Alli^ y The last rows of A and A! are equal. We can make both A and A! 
symmetric by multiplying by the last row the rows where the last column entry is X or Y. 
Thus the rows where A^i = X yield the same result as the rows where A!^ni — while we 


multiply the last row (which is equal in both matrices) by the rows where A^ j = A!l^ f. = Y. 
This means that A^ j^ = A'X = Z, so they do not get multiplied in their respective matrices. 
The effect of this is to complement the neighbourhood of node n, which comprises nodes j 
and k, and to induce the permutation {X Y) in the diagonal position of these rows. 

On the other hand, if we perform the local unitary operation ZK{DN)n, where Z = 

and ~ ^ J ^ largest set such that 

A^). = Z Wk E K, the effect is the same as described above (see |5] for a description of the 
effect of N). This implies that of the 6 parents, only at most 3 are local unitary inequivalent. 
□ 

We therefore obtain the following result. 

Lemma 25 Let A be such that, from Then A represents the mixed graph state 

5 5 

P = Y^ OTi, 5]] 9 = 1, Cj > 0, Vj. 


3=0 


3=0 


where pj = |0 oj)(0oj| + |0ij)(0ij|; 0 < j < 6, are the six density matrices obtained from 
the six parents as stabilized by A^y 


Proof. Follows immediately from lemma 23 □ 

Although only at most 3 of them are locally inequivalent, we include them all in the 
general sum p, since changing pj for a locally equivalent p) does not necessarily give local 
equivalence in p. 
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